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Materials containing hydrogen and high-Z ion species are used in a variety of high energy 
density physics (HED) applications, and it becomes increasingly important to investigate transport 
properties of these substances as they become fully-ionized multispecies plasmas. In particular, of 
significant practical interest, is the analysis of the plasma viscosity as it can determine stability and 
turbulence properties of the HED system. It has been pointed out by D. Ryutov [Phys. Plasmas 3, 4336 
(1996)] that the presence of a substantial amount of hydrogen in a multispecies plasma can significantly 
enhance its viscosity. The present work reports on a detailed analysis of this effect and the accurate 
calculation of the corresponding viscosity coefficient. 

Introduction-Materials, containing a considerable fraction of hydrogen or its isotopes 
(e.g., CH2, CD2, LiH) are utilized in a great variety of high energy density physics and inertial 
confinement fusion (ICF) applications. For instance, polystyrene (CH) or polyimide 
(C22H10N2O4) ablator shells are employed in the design of ICF capsules [1-2]. Deuterated 
polystyrene (CD) fibers have been considered for use as liners for wire array Z-pinches [3]. 
Local spherical blankets composed of the lithium-hydride (LiH) material can be potentially 
employed in ICF pulsed systems for tritium breeding and efficient (MHD) energy conversion [4- 
5]. In addition, CH2 plastic disk targets are extensively used in HED laboratory astrophysics 
experiments for laser-produced fast plasma flows [6-8]. It is therefore of significant practical 
interest to investigate the hydrodynamic transport properties of these materials as they become 
fully-ionized multispecies plasmas. 

The present work reports on the calculation of the viscosity coefficient, rj, which affects 
hydrodynamic stability and determines turbulence properties of the aforementioned HED 
systems. It has been previously pointed out by D. Ryutov [9] that the plasma viscosity can 
become enhanced by an order of magnitude due to the presence of hydrogen or a hydrogen 
isotope. Here, a quantitative analysis of this effect is performed, and an accurate expression for 
the multispecies plasma viscosity is derived. While the analysis is straightforward and is based 
on a well-known formalism [10-13], we are not aware of such calculation being done, and 
believe that the results of this work may prove useful in detailed studies of multispecies HED 
plasmas. 
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Qualitative analysis. -Consistent with the practical applications mentioned earlier, we 
consider highly-collisional dense plasmas where the mean free path of a hydrogen species, X HZ 
is small compared to the characteristic length scale for variations in the species flow velocities, 
L, 

A HZ /L«1. (1) 
In such highly-collisional multispecies plasmas, the hydrogen species (H) equilibrates its flow 
velocity to that of the heavier high-Z ions (Z) over a short collisional time scale, z c ~ X HZ /V TH , 

where V TH is the hydrogen thermal velocity. Considering here only frictional and viscous forces, 
the momentum fluid equations describing much slower relaxation of the common shear flow 
velocity V = V x (y)x are then given by 

d{M H n H V x ) V TH ,„ , d dV x 

" Hx = ~ -f~M H n H AV X + —7J H —± , (2) 
dt A HZ dy dy 

d(M 7 n 7 V ) V T7 d dV 

V z * = -r^M z n z AV x + ir f] z ^. O) 
dt A ZH dy dy 

Here, M, W, n, X, and n are the species mass, thermal velocity, number density, collisional mean 
free path, and viscosity coefficient, respectively, and AV* is the small difference between the 
hydrogen and high-Z species flow velocities. The viscosity coefficients in Eqs. (2) and (3), in 
general, include contributions from both the intra- and inter-species collisions. However, 
provided Z » 1 and assuming comparable hydrogen and high-Z species number densities, we 
neglect small contributions from weak H-H (relative to H-Z) and Z-H (relative to Z-Z) collisions, 
and estimate the collisional viscosity coefficients as rj H ~ rj HZ ~ M H V TH n H /l HZ and 

Vz ~ Vzz ~ M Z^TL n Z^ZZ • 

From Eqs. (2)-(3) it follows that slow relaxation of a shear in the common flow velocity 
occurs on a transport time-scale, T shear ~ (M Z n z /M H n H )(L/A HZ ) 2 t c » T c and is governed by 

d[{M H n H +M z n z )V x ]/dt = ^-Jj^- ~ ^-{M z V TZ n z A zz + M H V TH n H X HZ )^ . (4) 

dy dy dy dy 

It is also straightforward to show that the small difference between the hydrogen and high-Z 
species flow velocities is the order of AV x /v x ~ A 2 HZ /L 2 . Assuming equal temperatures of the 
hydrogen and high-Z ions, we readily obtain A HZ ~ Z 2 A ZZ , V TH ~ ^M Z /M H V TZ , and 
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Vhz n 



M, 



-^-^LZ 2 — tL. (5) 
Vzz n z V M z 

For instance, for a CD 2 plasma (n H /n z =2, M H /M z = 2/12 ,Z = 6), it follows that 
V hz I Vzz = 30 [9]. That is the hydrogen species provides the dominant contribution to the plasma 
viscosity, 

V « Vhz ■ (6) 
Note that the energy exchange between species occurs much faster, T E ~ (M z /M h )t c , than 

relaxation of the velocity shear, T c « T shear , and therefore the assumption of equal temperatures, 
T z = T H , used in deriving Eq. (5) is justified. 

While an accurate expression for the shear viscosity coefficient in a single-ion-species 
plasma is well-known [10], 

t. Jm t 5/2 

where lnA z denotes the Coulomb logarithm, the goal of the present work is to derive an 
equivalent accurate expression for the viscous coefficient corresponding to the H-Z collisions, 
i.e., 7] HZ . 

Quantitative analysis.-lt is well-known that the viscosity of a monatomic gas or 
unmagnetized plasma is described by a single shear-viscosity coefficient, 7], with the viscous 

stress tensor given by a tj = -rfi^V^dXj +dV j jdx i -2/3^-V • v) [10-11]. From the standard 

framework for determining transport coefficients [10-12], it follows that in order to evaluate the 
viscosity coefficient, rj HZ , one needs to find a first-order correction in A, HZ /L « 1 , i.e., Fi, to a 
zero-order solution of the kinetic equation given by a local Maxwellian distribution with a shear 
flow velocity 



71' V' 



TH 



[u x -vAyf + u 2 y+ u^ 



(8) 



Here, r and u are the particle's coordinate and velocity, V TH = y/2T H /M H is the hydrogen 
thermal velocity, and, without loss of generality, we adopted Y = V x (y)x. Introducing 
u x =u x —V x (y), one obtains the following equation for the first-order correction F\ 
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U X U V dV r i 

1 -hT-^r- F u=v B z<fr l l (9) 

V m dy 

where, v HZ (u) = {l7m z Z 2 e A / M 2 H u 3 )ln A H , and C[f] denotes a normalized collision operator for 
H-Z collisions. In what follows we consider a succession of increasingly detailed collision 
operator options. First, we assume a small mass ratio M H /M Z «1, and adopt the Lorentz 
collision operator to obtain an exact solution to Eq. (9). A first-order correction in M H /M z « 1 
is then analytically recovered by assuming a Maxwellain distribution of the high-Z ions and 
considering hydrogen collisions with this Maxwellian background. Finally, the same collision 
model is used to obtain numerical solutions to Eq. (9) for the case of an arbitrary mass ratio. 

Case (a): Lorentz model. -Assuming, for simplicity, M H «M Z we adopt the Lorentz 
collision model, 

d . „dE 1 d 2 F, 



r n Id. „dF, 
C[Fj = - sin# — - 



+ -^^, (10) 



sin#30" '"" dd sin 2 0 d(p 2 

where, (u,0,<p) are the spherical coordinates corresponding to the coordinate system [u x ,u y ,u,). 

In general, the collision operator in Eq. (10) should also include terms proportional to the relative 
velocity between the H and Z species (see Refs. [10]-[12]). However, these terms are of the 

second order in A HZ /L (recall that AV x /v x ~ A 2 HZ /L, 2 ), and therefore are neglected in Eq. (10). 
Noting that 

u x u y /u 2 = cosff7sinflcosfl = 2 "'^ — ~. (H) 

where Y /m (^,^?)are the spherical harmonic functions [14] that satisfy c(y, m )=-l(l + l)Y, m , it 
follows that 

F fS^ Y^-F,, u 2 F M dV x 

~ V77 2 77~2 r^^T - • ( 12 ) 

V 15 6 V 2 H v HZ {u) dy 



The viscous momentum flux is given by 

— M f 77 i, 77 rl 3 n — —n 

dy ' 



= v yx = M H J u x u y F x d u = -J] HZ — i , (13) 



and we readily obtain 



rj HZ =M H §dn(Y 2 ,_ l -Yj 2 ]du^x^^ (14) 

4 j V TH V HZ \U ) 



where dQ. - sin 6d0d(p . Making use of the properties of spherical harmonics, i.e., 
Y ; ,_ m =(-irY,; and f dfl^ = S lf S^ , it follows that \dd^ 2X -Y 2 J =2. 
Performing the remaining integration along the u-axis, we obtain 

32V2 n H ^M~ H TT 



Vhz = 



or, equivalently, 



I5x 3 ' 2 n z zVlnA H 



,1/2 



(15) 



In A, 



(16) 



where A H is the atomic mass of the hydrogen isotope. Recalling Braginski's result for r] zz [Eq. 
(7)], a more accurate form of the estimate in Eq. (5) is given by 



HZ 



= 1.33^-Z 



M, 



V 



zz 



M, 



(17) 



where lnA z = lnA w is assumed for simplicity. It is also straightforward to extend the result in 
Eq. (15) for the case where several high-Z ion species are present 

32V2 n H JM^T'J 2 



Vhz = 



\5n 312 £Z> Z . e 4 lnA„ 



(18) 



Case (b): Collisions with a Maxwellian background, M H /M z « 1 .-The expressions for 
the viscosity coefficient in Eqs. (15)-(18) are only valid in the limit M H «M Z , and do not 
include any effects associated with a finite mass ratio. A first-order correction that includes the 
effects appearing in order 0(M H /M z ) can be readily derived as follows. Assuming that the 
charge state of the Z-species ions is sufficiently high, and thus Z-Z collisions are sufficiently 
strong, we neglect deviations of the high-Z species distribution from a Maxwellian. The collision 
operator in Eq. (10) describing H-Z collisions can then be generalized to include effects of an 
arbitrary mass ratio as follows [12] 



C[F x ] = v D 



sin 6? 



3F, 



+ 



1 d 2 K 



sin 



odd de %m 2 ed(p 2 



+ u- 



du 



M H r 1 dF, 



K M H +M Z 



du 



• (19) 



Here, V D = Y{g)-G{g), V s = 2[{M h + M z )/T z ]g(^)/u and v n =4G{^)/u 3 are the collision 
coefficients, % = uj ^2T Z /M ~ is the normalized particle speed, and the dimensionless functions 
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y(#) and G{£) are specified by Y{£)= )jV' 2 <fr and G{£) = [y{£)- xY '{£)]/ 2% 2 . 

Although the collision model in Eq. (19) is valid for an arbitrary value of M H / M z , to progress 
analytically here we assume M H /M z « 1 . The operator in Eq. (19) then takes on the following 
form: 



V M 



2u M 



z J 



-sin# 



dE 



+ ■ 



1 d 2 F^ 



v sin 0 30 dd sin 2 <9 d(p 2 , 



+ 2 — —u — 

M z du 



^ T 7 1 dE 

F 1 +— L 

M H u du 



■ (20) 



Representing the solution to Eq. (9) as Fj = F° + (M H /M z )f/ 1 ', where Fj° is given by Eq. (12), 
the correction f/ 1 ' is determined from 



-sin6> — + - 



»2 A 



sin 6d6 d0 sm 2 6d(p 2 



A 



p(l) V TH M H pQ 

2u 2 M, ' 



+ 2u- 



J 



du 



M + 



M „ u du 



0. (21) 



Recalling that equilibration of species temperatures occurs much faster than relaxation of a 



common velocity shear, we take T z =T H = T , and after some straightforward algebra obtain 



Vhz = 



32V2 n H y[M^T 5/2 



I5x 3/2 n 7 Z 2 e 4 \nA, 



1 11 M H 



12 M 



(22) 



z J 



l Z ^ " lll -"-H 

In the case several high-Z ion species are present, the viscosity in Eq. (22) takes the following 
form 



Ihz = 



32V2 



n, 



5/2 



15^ 3/2 £Z> ZJ e A \nh H 



ll^M 



12 y« 5 



z/ 



(23) 



Case (c): Collisions with a Maxwellian background, M H /M Z ~1.-The expressions in 
Eqs. (22)-(23) omit higher-order corrections in M H /M Z [i.e., C>(m 2 /M 2 ), o(m* /M*), ...]. 
However, these terms may be important for a detailed analysis of the lithium-hydride HED 
plasmas that will appear in an ICF reactor if a lithium-hydride blanket is used [4-5]. In order to 
address this practical issue, which is also of general theoretical interest, we now evaluate 
accurate values of the H-Z viscosity coefficient for an arbitrary M H /M z mass ratio, by making 

use of the variational principle [13]. Introducing /, =FjF M , where F l is an arbitrary trial 
function, and noting that the collision operator in Eq. (19) is self adjoint with 

I d V,v HZ (u)c[f; ] = J d \F iVhz {u)c[l ] , (24) 
it follows that Eq. (9) is the Euler equation for the variational principle 
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3S = S{S l -2S 2 ) = 0, (25) 
where the functionals Si and S2 are defined by 

S x =\d'xxf x v HZ {u)c[FX (26) 

r o - 2u u dV 
S 2 = \d 3 uf 1 -f 2 -^-±F M . (27) 



It is important to note that the extremal value of the variational quantity, S 6 " 1 , satisfies 

(28) 



gext ^ext 



' 2 BV^ 
K M H V T 2 H dy j 



where S ex ' and S e 2 xt corresponds to the functionals in Eq. (25) and Eq. (27) evaluated with the 
exact solution to Eq. (9). The fact that the viscosity coefficient is directly related to the extremal 
value of the variational quantity provides very good accuracy of this asymptotic method. 

We now proceed with the variational analysis by making use of the following trial 
function F x [10] 

F y (u) = [Y 2M {6, <p) - Y v {6, <p% 2 X a,.Lf /2) ( tf 2 /V T 2 H ) , (29) 

i=l 

where l[ 5/2) {x) are Sonine polynomials, L^ /2) {x) = l, l[ 5/2) {x) = 7/2- x, 
L* 2 5/2 ' (x) = 63/8 - 9x/ 2 + x 2 /2 , N determines the number of terms retained the expansion, and the 
coefficients a, are found from a linear system of equations, dS/da i - 0 . Due to orthogonally of 
the Sonine polynomials, the viscosity coefficient, rj HZ , is determined only by the value of ao. 
While an exact value of ?] HZ would be recovered in the limit N — > °° , it is found that sufficient 
accuracy can be obtained by retaining only a few first terms in the expansion in Eq. (29). The 
results for the viscosity coefficient can be expressed as 



7m= ^^^l!L (l _ B ), (30 ) 
15^ 3/2 n z Z 2 e 4 lnA H 

where the value of the coefficient B is given in Table 1 for various choices of M H /M z . It is 
important to emphasize here that the result in Eq. (30) was obtained by making use of the 
collision operator in Eq. (19), i.e., neglecting deviations of the high-Z species distribution from a 
Maxwellian. Also, equal temperatures of hydrogen and high-Z species were assumed. 
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Table 1. Numerical values of the coefficient B [in Eq. (30)] obtained with N=4. 



MzJM H 


1 


2 


3 


4 


5 


6 


7 


8 


9 


»1 


B 


0.374 


0.283 


0.220 


0.178 


0.150 


0.129 


0.113 


0.101 


0.091 


11 M H 

12 M z 



Note that unlike the results in Eqs. (15) and (22), the result in Eq. (30) cannot be 
extended straightforwardly to the more general case of several high-Z ion species. Each new 
combination of high-Z species requires separate use of the variational principle to determine the 
corresponding viscosity coefficient. However, we emphasize that the first-order accurate 
estimate of the viscosity coefficient for the case where several heavy high-Z ion species with 
M H /M Zj « 1 are present is given in Eq. (23), or equivalently 



Vhz( 



i -i 
s cm g 



)= 



4.27 xlO" 5 n H Af{T[eV]) 



5/2 



In A, 



2» 



Zj 



ll^M 



Zj 



(31) 



where A H is the atomic mass the hydrogen isotope. 

In conclusion, as an illustrative example, we present numerical values of the viscosity 
coefficients for materials that are often considered for HED and ICF applications (see Table 2). 



Table 2. Viscosity coefficients for several choices of a hydrogen-containing materiaf , 

axlO" 6 



T] HZ (s~ 1 cm~ 1 g)=- 



In A, 



■(W 



5/2 



Material 


CH 


CH 2 


CD 


CD 2 


LiH 


C22H10N2O4 


a 


1.10 


2.20 


1.56 


3.12 


4.20 


0.346 



a The coefficient a is obtained by using Eq. (31) for the case of polyimide (C22H10N2O4), and by 
using the variational method for all other materials. 
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